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Abstract. After reviewing some basic relevant properties of stationary stochastic processes (SSP), 
we discuss the properties of the so-called Harrison-Zeldovich like spectra of mass density pertur- 
■ bations. These correlations are a fundamental feature of all current standard cosmological models. 

Examining them in real space we note they imply a sub-poissonian normalised variance in spheres 
Om(R) ~ R 4 lnR. In particular this latter behaviour is at the limit of the most rapid decay (~ R~ 4 ) 
of this quantity possible for any stochastic distribution (continuous or discrete). In a simple clas- 
sification of all SSP into three categories, we highlight with the name "super-homogeneous" the 
properties of the class to which models like this, with P(0) = 0, belong. In statistical physics lan- 
guage they are well described as lattice or glass-like. We illustrate their properties through two 
simple examples: (i) the "shuffled" lattice and the One Component Plasma at thermal equilibrium. 
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^ INTRODUCTION 

In standard theories of structure formation in cosmology the density field in the early 
Universe is described as a perfectly uniform and isotropic matter distribution, with su- 
perimposed tiny fluctuations characterized by some particular correlation properties (e.g. 
[[p). These fluctuations are believed to be the initial seeds from which, through a com- 
plex dynamical evolution, galaxies and galaxy structures have emerged. In particular, 
in all the standard models, the initial fluctuations are taken to be Gaussian and with a 
power spectrum (PS) P(k) satisfying the so-called Harrison-Zeldovich (HZ) condition 
of being proportional to k at small k @]. 

The present paper has two main porpuses [|3|]. Firstly, to clarify the real space statisti- 
cal properties of the mass density fluctuations common to all the standard cosmological 
models (which have been almost completely overlooked in the literature on the sub- 
ject). And secondly, through this discussion, to relate and compare these models of the 
primordial Universe to correlated systems encountered in statistical physics. 

In particular we find that all these standard cosmological models are characterized 
by a "superhomogeneous" (or superuniform) matter distribution. This means that mass 
fluctuations over sufficiently large spatial scales are "sub-poissonian", i.e. they increase 
with the spatial scale more slowly than in a random poissonian matter distribution. 
Moreover, it is shown that the scaling of mass fluctuations satisfying the HZ condition 
approaches the slowest possible for any stochastic matter distribution. We will see that 
in the context of usual statistical physics this kind of behavior is recovered in the case of 



lattice-like or glass-like particle distributions or in the so called One Component Plasma 
(OCP) [Q] at thermal equilibrium. In this context we present some simple recipe to build 
particle distributions satisfying the HZ condition at small k. This last point can be useful 
first of all in the context of Af-body simulations for the study of structure formation from 
primordial fluctuations through a gravitational dynamics 



BASIC CONCEPTS OF CORRELATION ANALYSIS 

Inhomogeneities of the mass density field in cosmology are described using the gen- 
eral framework of stationary stochastic processes (hereafter SSP). Let us consider in 
general the description of a continuous or a discrete homogeneous mass distribution 
p(r) in terms of such a process. A stochastic process is completely characterized by its 
"probability density functional" £?\p{j)\ which gives the probability that the result of 
the stochastic process is the density field p(r) (e.g. see Gaussian functional distributions 
[^p). For a discrete mass distribution, the space (e.g. the infinite three dimensional space) 
is divided into sufficiently small cells and the stochastic process consists in occupying 
or not any cell with a point-particle following certain correlated probabilities, and p (r) 
can be written in general as: 

oo 

P(r) = £S(r-r ; .), 
i=l 

where r- is the position vector of the particle i of the distribution . 

The word "stationarity" refers in the present context to the spatial stationarity of the 
process, and means that the functional &[p(r)] is invariant under spatial translation. We 
suppose also that the distribution is statistically isotropic (invariance of J^[p(r)] under 
spatial rotation), and has a well defined positive average value: 

(P(?)>=P >° 5 

where (...) is the ensemble average over all the possible realizations of the stochastic 
process, i.e. the average over the functional ^[p(r)]. Moreover it is usually assumed 
that 2?\p(j)\ is ergodic. This means that spatial averages in single infinite realization 
coincide with the ensemble averages. 

The quantity (p(r ] )p(r 2 )...p(r l )') is called the complete /-point correlation function. 
In the discrete case (p(r 1 )p(r 2 )...p(r l ))dV 1 ,dV 2 , ■■■,dV l gives the a priori probabil- 
ity of finding / particles, in a single realization, placed in the infinitesimal volumes 
dV l ,dV 2 ,---,dV l respectively around r v r 2 ,...,r l . In both cases the statistical stationar- 
ity and isotropy (hereafter SSI) imply that the /-point correlation functions, for any /, 
depend only on the scalar relative distances among the / points 

Let us analyze in further detail the auto-correlation properties of these systems. As 
aforementioned, in our hypothesis, (^p(r ] )p(r 2 )^ depends only on r 12 = \ f\ — f 2 \. 

The reduced two-point correlation function | (r) is defined by: 

{p{r x )p{f 2 ))=pl 



1 +f(r 12 ) 



The correlation function f(r) is one way to measure the "persitence of memory" of 
spatial variations in the mass density [£p. For a discrete distribution of particles this 
means that p <^(r)J 3 r measures the excess of probability of finding a particle in a 
small volume d 3 r at a distance r from another fixed particle, with respect to the a 
priori probability p Q d 3 r. Obviously, since p must be the average density in any single 
realization in the infinite volume limit, f (r) satisfies the condition of vanishing in the 
limit r — > °°. 



The Poisson particle distribution 

Let us consider the paradigm of a stochastic homogeneous point-mass distribution: 
the Poisson case (see [0]). It can be defined as follows: let us partition the space in cubic 
cells of volume AV, and then occupy each cell independently of the others, with a point- 
particle with a probability n Q AV, where n Q > and AV <C 1 . It is simple to show that 
in the limit AV — > 0, one can write: 

Po = n o (!) 
I(r)= 8 -^ (2) 

where 8(r) is the usual Dirac delta function. Equation @ is a direct consequence of 
the fact that there is no correlation between different spatial points in the definition 
of the stochastic process. That is, the reduced correlation function ^ has only the so 
called diagonal part. It is simply shown that this diagonal part is present in the reduced 
two-point correlation functions of any statistically homogeneous discrete distribution of 
particles with correlations. Therefore in general [^p for a SSI distribution of particles the 
reduced correlation function can be written as 

f, s 8(r) , 
Po 

where t, is the non-diagonal parts which is meaningful only for r > 0. 



THE MASS VARIANCE IN A SPHERE 

In this section we consider the amplitude of the mass fluctuations in a generic sphere of 
radius R with respect to the average mass. First let M(R) = J c ^p(r)d^r be the mass 

(for a discrete distribution the number of particles) inside the sphere C(R) of radius R 
and then volume \\C(R) \\ = ^j-R 3 - The normalised mass variance is defined as 

2 (M(R) 2 } - (M(R)) 2 



It is simple to show that o^(R) in a stationary stochastic mass density can be rewritten 



as: 



This formula will be useful in the following for a complete classification of the mass 
density fields with respect to the large scale behavior of mass fluctuations. 
Note that since | (r) — > for r — > °°, then 

Jim oliR) = 0. 

This is nothing but the condition of existence of a well defined average density p . 



A classification of mass fluctuations 

If we apply Eq. @ to the case of a Poisson point process as above defined, we find 

2 1 _ 1 

° m{r) ~ = WW) ' 

that is (AM 2 (_/?)) ~ \\C(R)\\. The same widespread behavior is found in systems char- 
acterized by mainly positive short range correlations (e.g. a perfect gas at high tempera- 
ture). 

We can characterize an arbitrary mass distribution with respect to the behavior of mass 
fluctuations as follows. Let us suppose o^(R) ~ R a at large R, then: 

• If < a < 3 the system has critical fluctuations, typical of the order parameter of a 
thermodynamical system at the critical point of a second order phase transition. As 
we show below, this happens in the case of mainly positive long range correlations 
such that / d 3 rt; (r) = +°°; 

• If a = 3 we have a substantially poissonian system, characterized by mainly posi- 
tive short range correlations, such that / cPrE, (r) = c where < c < +°°; 

• If a > 3 the system has sub-poissonian fluctuations typical of lattice or glass-like 
systems with an almost ordered arrangement of mass perturbations (of particles in 
the discrete case), characterized by a balance between positive and negative cor- 
relations, such that / d 3 r<!; (r) = 0. For this reason we call these systems superho- 
mogeneous. We show also below that a < 4 (a < d + 1 in d dimensions) for any 
stochastic mass distribution. As shown below all the standard cosmological models 
for primordial mass fluctuations belong to this class with a = 4 (for a difference 
between a Poisson and a superhomogeneous distribution see Fig. [[]). 




THE POWER SPECTRUM P(K) 

Much more used in cosmology than | (r) is the equivalent £-space quantity, the power 
spectrum P(k) which is defined as 



P(k) = lim 



S P (k)\' 



where 8 p (k) = J v d re r (p (r) — p ) / p is the Fourier integral in the volume V of the 
normalized fluctuation field (p(r) — p )/p . In a SSI matter distribution this depends 

only on k = \k\ and can be written as 

P(k) ee P(k) = I d Z re- llrr l(r) = An T Jrr 2 ^^|(r) . (4) 
J 'Jo kr 

It follows from its definition that P(k) > (this property and Eq. (Q) constitute the so- 
called Khincin theorem of SSP's). In particular for a Poisson particle distribution one 
has P(k) = p l . At last note the important relation P(0) = / J 3 r| (r) for any SSP. 



The cosmological power spectrum, and the HZ condition 

In current cosmological models it is generally assumed that the power spectrum 
P(k) ~ k at small k. This is the famous Harrison-Zeldovich (HZ) condition. It is believed 
to describe the "primordial" fluctuations at very early times, the putative remnants of a 
period of "inflation" prior to the ordinary Big Bang phase [jlj, [TT]]. The HZ behaviour is 
appropriately cut-off at a short distance scale (i.e. for k larger than some cut-off k c ). The 
type of the cut-off depends on the single cosmological model (e.g. Cold Dark Matter, 
Hot Dark Matter, etc.). The reason for the HZ condition is tied to considerations about 
the consistency of the Friedmann-Robertson- Walker (FRW) metrics under perturbations 
(see [[]], HO)]): any other spectrum will give mass fluctuations which dominate over the 
homogeneous background either at some time in the future or past. 



One can see [|3|] that in all the standard cosmological models the corresponding f (r) 
is positive at small r and has a negative tail ~ — r~ 4 at large r with only one zero. 



The power spectrum vs. the mass variance 

Quite generally, by studying Eq. (Q), it is not difficult to show [BR that for a spec- 
trum P(k) ~ k n for — > (note that the Poisson-like system are given by n = 0) and 
appropriately cut-off at large fc, one has at large R 

f X l R 3+n ■ l f n< \ 

o 2 M {R)~{ log(R)/R 4 if n = \ (5) 
[ 1 /R A if n > 1 

In terms of the non-normalized quantity (AM 2 (_/?)) ~ o^j(R)R 6 , Eq. (|5|) says that for 
n>0(i.e. /d 3 r|(r)=0) we have a scaling behavior of (AM 2 ) as a function of 7? slower 

than for Poisson fluctuations, which correspond to n = (i.e. / rf 3 r| (r) = c > 0), with 
the limiting behavior for n > 1 corresponding to quadratic mass fluctuations which are 
proportional to the surface area of the sphere. These systems are thus characterised by 
surface fluctuations, ordered (or homogeneous enough, one could say) to give this very 
particular behaviour. The case n < (i.e. / J 3 r^ (r) = +<*>) is typical of critical phenom- 
ena in which long range mainly positive correlations determine a "super-poissonian" 
behavior of quadratic fluctuations. 

Equation (Q) shows clearly that cosmological models, characterized by the HZ condi- 
tion (n = 1), are not only superhomogeneous, but are at the limit of the slowest possible 
behavior of quadratic mass fluctuations with the spatial scale. 

TWO "SUPERHOMOGENEOUS" EXAMPLES 

In this section we present two examples of superhomogeneous point-particle distribu- 
tion: the shuffled lattice [g] and the One Component Plasma (OCP) [[§. 

1) In the case of particles placed on the sites of a regular cubic lattice (in any 
dimension d), one has [|3|] p = aT d where a is the lattice spacing and: 

where the sum is extended to all the vector of the reciprocal lattice but the origin 
H = 0. Therefore we can say that around k = Pj(k) ~ k + °°, and then it is the most 
superhomogenehous particle distribution with o^R) ~ R~( d+l \ 

At this point let us introduce a stochastic displacement field, in which each particle 
is displaced from its lattice site R of a random and statistically isotropic vector u~, each 
particle independently of the others. One can show [ |TT| ] that, if the probability density of 
the displacements p(u) has a finite variance, then the resulting particle distribution has 



a power spectrum P sl (k) ~ k 2 at small k. If instead, the variance is infinite, and p(u) has 
a tail going as uT a (with J < a < d + 2 for normalizability, but with infinite variance), 
then P i(k) ~ k a ~ d at small k. Therefore, choosing appropriately a we can obtain any 
kind of superhomogeneous particle distribution with the exponent of P sl (k) at small k 
ranging from to 2. In particular for a = d + 1 (a = 4 in three dimension) we obtain the 
HZ condition. 

2) Let us now present briefly another model taken from equilibrium statistical physics 
displaying superhomogeneous mass fluctuations: the OCR It is a system of interact- 
ing point particles carrying a unit positive charge and mass repulsing each other via a 
Coulomb potential V(r) = l/r, in a continuous uniform background negatively charged 
giving overall charge neutrality [f?p. One can show quite generally that, at thermody- 
namic equilibrium the system has only one phase, and at small k the power spectrum of 
the particle distribution satisfies: 

a 

P(k) 



4xp$P ' 

where p is the average density of particles and /3 is the usual inverse temperature. It is 
very interesting to note that [O] if, instead of having the Coulomb potential V(r) ~ r _1 , 
one has a potential V(r) ~ r one would obtain 

P(k) 



2n 2 n 2 [5 



that is the HZ condition. These results can be very important if applied to A'-body 
simulations of gravitational dynamics to study structure formation starting from the 
primordial mass fluctuations In fact, in this context a great importance is given to the 
accurate preparation of the initial conditions for the simulation. For instance, this can be 
set through a modified OCP with V(r) ~ r~ 2 at large scales and an appropriate modified 
behavior at smaller scales [[T2p . 

Both 1) and 2) provide particle distributions similar to the right side picture of Fig. []]. 
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